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Abstract 


The  use  of  plasma  ignition  for  high-performance  guns  has  been  pursued  at  several  levels.  As 
part  of  the  program  at  the  U.S.  Army  Research  Laboratory  (ARL),  we  are  attempting  to  define 
the  role  of  the  broadband  electromagnetic  radiation  that  is  a  key  energy  component  of  the  plasma 
discharge.  In  particular,  this  effort  is  developing  a  radiation  transfer  model  that  will  be 
incorporated  into  the  computational  fluid  dynamics  (CFD)  model  of  plasma  ignition.  This  effort 
will  allow  more  physically  realistic  characterization  of  the  energy  transfer  in  the  plasma  and  its 
interaction  with  the  propellant  in  the  models. 

This  report  summarizes  the  physical  background  of  the  radiation  transport,  and  the  validation  of 
the  different  submodules  by  use  of  published  data  will  be  presented.  Current  efforts  include  the 
coupling  of  the  radiation  model  to  the  ARL  NSRG2  code  for  the  plasma-propellant  interaction 
simulation.  Furthermore,  a  series  of  laboratory  experiments  will  be  performed  to  validate  the 
CFD  code  with  the  incorporated  radiation  transport  submodel. 
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1.  Introduction 


The  plasma  ignition  of  high-loading  densities  has  been  investigated  for  many  years  and  in 
various  gun  calibers.  The  advantages  of  these  adjustable  igniters  (e.g.,  a  temperature 
compensation  effect)  have  been  demonstrated  in  many  firings,  and  the  interaction  of  the  plasma 
with  the  propellant  has  been  the  focus  of  many  open-air  and  closed- vessel  experiments.  Yet,  not 
all  of  the  physical  and  chemical  effects  of  the  plasma-propellant  interaction  (PPI)  are  understood 
in  adequate  detail. 

Many  experimental  investigations  have  indicated  that  the  energy  transfer  by  plasma  radiation  is 
an  important  effect  for  a  fast  ignition  of  the  propellant  as  well  as  an  enhanced  mass  conversion 
[1-4,  5],  Whether  or  not  it  dominates  all  other  energy  transport  mechanisms  from  the  plasma 
source  to  the  propellant  will  depend  strongly  on  the  ignition  concept  (e.g.,  Plasma-Jet,  Current 
Injection)  and  its  parameters  (like  geometry,  pulse  duration,  pressure  levels,  etc.).  Independent 
of  the  PPI  effects,  radiation  is  the  dominant  effect  of  energy  transport  within  the  hot  plasma.  It 
has  been  known  for  many  years  that  the  ohmic  heating  of  the  plasma  (by  the  electrical  energy 
release)  is  mainly  balanced  by  radiation.  Therefore,  radiation  transport  has  to  be  taken  into 
account  when  solving  the  energy  balance  equation  in  any  electrothermal  chemical  (ETC) 
computational  fluid  dynamics  (CFD)  calculations  and  when  studying  and  simulating  PPI  effects. 

In  this  report,  the  status  of  a  new  radiation  transport  model  will  be  described,  which  will  be 
coupled  to  the  ARL  CFD  code  NSRG2  for  a  more  physically  realistic  description  of  plasma  in  it. 
The  coupled  codes  then  will  be  used  to  study  the  PPI  and,  finally,  to  improve  the  overall  interior 
ballistic  simulation  of  ETC  guns.  In  the  hot  plasma,  the  radiation  by  atoms  and  ions  (the 
dissociation  and  ionization  products  of  the  capillary  material  or  eventually  of  the  propellant) 
dominates  over  molecules,  and  up  to  now  we  have  focused  on  this  hot  radiation  source.  But  in 
the  next  step,  the  absorption  by  the  most  important  molecules  (in  the  cold  boundary  layer  of  the 
plasma  and  the  surrounding  cold  propellant  gases)  will  be  implemented  as  well.  An  adequate 
description  of  the  radiation  transport  code,  with  all  of  its  physical  background,  would  go  far 
beyond  the  scope  of  this  report;  only  a  summary  of  the  basic  theories  with  a  focus  on  some 
important  aspects  is  presented  next. 


2.  Physical  Background 


In  this  section,  the  physics  of  plasma  radiation  is  summarized;  the  reader  is  referred  to  some 
earlier  work  [5,  6],  the  final  report  of  this  Short-Term  Analysis  Service  (STAS)  contract  [7],  and 
the  references  therein  for  further  details.  The  computation  of  the  radiation  transport  equation 
(RTE)  depends  on  frequency-dependent  absorption  coefficients  kv  and  emission  coefficients 
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which  describe  the  interaction  of  the  radiation  field  with  matter  (i.e.,  the  atoms,  ions,  and 
molecules).  The  relation  between  absorption  and  emission  are  given  under  the  assumption  of  a 
local  thermodynamic  equilibrium  (LTE)  by  Kirchhoff  s  law: 

sv=kvIqv,  (1) 

where  7®  is  the  intensity  of  a  black  body  radiator  at  frequency  v  and  temperature  T.  Therefore,  it 
is  sufficient  to  calculate  either  the  emission  or,  more  commonly,  the  absorption  coefficients,  if 
LTE  is  assumed.  These  absorption  coefficients  depend  on  cross  sections,  CTijk,  for  line  and 
continuum  absorption  and  on  the  number  densities,  Nyk,  of  species  k  in  ionization  stage  j  and 
excitation  level  i.  Some  aspects  in  the  theoretical  calculation  of  these  quantities  will  be 
discussed  in  the  following  subsections. 


2,1  Equation  of  State  (EOS) 


An  essential  difficulty  in  obtaining  equations  of  state  for  plasmas  arises  from  the  fact  that  the 
Coulomb  potential  of  the  ionized  species  give  rise  to  an  infinite  number  of  energy  levels.  The 
high  levels  have  large  radii,  which  eventually  become  larger  than  the  mean  interparticle  distances 
in  the  plasma.  Furthermore,  the  required  partition  functions  then  become  infinite.  For  these 
reasons,  some  cut-off  procedures  have  been  introduced  in  the  past.  The  most  common  ones  are 
based  on  the  Debye  or  Debye-Hueckel  Theory,  which  leads  to  a  reduction  of  the  ionization 
potential.  At  higher  pressures,  when  the  Debye  Theory  is  not  valid  anymore,  often  the  Ion- 
Sphere  Model  is  introduced — all  levels  with  a  radius  larger  than  the  mean  interparticle  distance 
cannot  exist  in  this  theory.  Theories  which  assume  a  fixed  cut  off  at  any  exited  quantum  level 
often  lead  to  useful  results  when  calculating  thermodynamic  properties  of  hot  gases,  but  they 
lead  to  some  unphysical  behavior  when  dealing  with  optical  quantities,  i.e.,  the  sharp  edges  of 
photo-ionization  cross  sections  at  the  threshold  of  the  photo-ionization  energy,  which  never  have 
been  observed  in  experimental  measurements.  Within  the  framework  of  the  OPACITY  Project 
[8],  an  international  collaboration  concerned  with  the  calculation  of  opacities  of  stellar 
envelopes,  Hummer  and  Mihalas  [9]  introduced  an  occupation  probability  formalism  for  the 
calculation  of  the  internal  partition  function: 
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Here,  gyk  is  the  statistical  weight,  wjk  is  the  occupation  probability,  Eyk  is  the  energy  of  the  level, 
k/i  is  the  Boltzmann  constant,  and  T  is  the  temperature.  The  Wyk  go  to  zero  for  higher  states, 
which  ensures  that  the  sum  in  equation  2  is  convergent.  Interaction  with  charged  and  neutral 
surrounding  particles  is  taken  into  account.  The  perturbations  by  neutral  species  are  based  on  an 
excluded  volume  (hard-sphere)  treatment  (first  term  on  the  right  of  equation  3)  and  perturbations 
by  charged  particles  are  calculated  from  a  fit  to  a  quantum  mechanical  Stark  ionization  theory 
(using  a  first-order  perturbation  theory  and  a  Holtzmark  distribution  for  the  ion  microfield 
[second  term  of  equation  3]): 
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The  index  v  runs  over  all  neutral  species  only;  r\  v  is  the  ground  state  radius  of  the  neutral 
perturber,  and  is  the  radius  associated  with  state  i  of  species  k  in  ionization  stage  j  of  the 
perturbed  atom.  Zjk  is  the  net  charge  {Zjk  =  0  for  neutrals),  Xijk  is  the  ionization  energy,  and  kijk  is 
a  quantum  mechanical  correction  and  the  sum  in  term  two  runs  over  all  charged  particles 
excluding  electron.  By  use  of  this  occupation  probability  formalism,  the  number  densities  in 
each  excitation  state  are  then  given  by  a  modified  Boltzmann  equation: 

Myk  =  N  jkgijkV  ^  'kBT  /  Z’jk  •  (4) 


A  comparison  of  the  simulated  emission  spectrum  of  the  hydrogen  B  aimer  series  to  experimental 
results  by  Wiese  et  al.  [10]  will  be  shown  later.  It  can  be  seen  that  the  occupation  formalism 
leads  to  dissolution  of  the  higher  lines  of  the  spectral  series,  resulting  in  a  smooth  transition  to 
the  continuum  and  an  excellent  agreement  with  the  experimental  measurements. 

For  the  calculation  of  the  ionization  equilibrium  number  densities,  an  extended  version  of  the 
NASA-Lewis  code  [11,  12]  is  used,  which  has  been  modified  to  simulate  hot  gases  and  plasmas 
[13,  14]  (extension  to  high  temperatures,  including  a  collision  cross-section  database  for 
calculation  of  transport  coefficients  and  using  a  virial-type  EOS  for  high-pressure  conditions).  The 
code  is  based  on  the  free-energy  minimization  procedure  in  the  “chemical  picture,”  assuming  an 
LTE.  In  the  future,  the  particle  number  densities  might  be  estimated  with  Nusca’s  NSRG2  code 
[15],  including  Anderson’s  rate  equation  submodule  [16].  Figure  1  shows  calculated  number 
densities  for  polyethylene  plasmas  at  temperature  and  pressure  conditions  that  typically  occur  in 
current  ARL  plasma  jet  experiments  [4]. 


2.2  Absorption  Coefficients 

The  cross  section  for  line  absorption  (bound-bound)  is  given  by 


me 


(5) 


with  the  charge  of  an  electron  e,  its  mass  me,  and  the  velocity  of  light  c;  f.r  is  the  oscillator 
strength  for  a  transition  i  — >  i' ,  and  IU'  is  the  profile  of  the  line. 

Line  profiles  for  atoms  and  ions  in  a  plasma  are  determined  by  electron  impacts  and  by  ion 
microfields.  Griem  [17]  has  made  extensive  calculations  for  these  processes,  and  Konjevic  et  al. 
[18]  have  reviewed  experimental  data.  It  was  found  that  the  experimental  data  are  in  good 
agreement  with  calculations  made  using  the  methods  of  Griem.  In  the  semi-empirical  electron 
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(a) 


(b) 


Figure  1.  Calculated  LTE  number  densities  of  a  polyethylene  (C2FI4)  plasma  at  (a)  1  MPa  and  (b)  10  MPa. 


impact  approximation  by  Griem  [17,  19]  and  Baranger  [20],  the  full  width  of  the  line  at  half 
maximum  (FWHM)  is  given  by 
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Here,  Rv  is  the  Rydberg  energy,  h  is  Planck’s  constant,  a0  is  Bohr’s  radius,  and  the  matrix 
elements  can  be  expressed  as 

\<i'\r\i>\2=S(J\J)/gi,  (7) 


with  the  line  strength  S  for  a  transition  g,  in  equation  7  is  the  statistical  weight  of  the 

ground  level,  and  the  sums  running  over  all  levels  with  an  allowed  transition  to  the  state  i  (/). 
The  g  values  in  equation  6  are  the  so-called  effective  Gaunt  factors  [  1 7]  for  the  initial  and  final 

states 
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The  profile  of  the  line  in  the  case  of  electron  impact  broadening  has  a  Lorentzian  shape 
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and  when  Doppler  broadening  is  also  taken  into  account,  the  convolution  of  the  two  profiles  is 
assumed  to  have  a  Voigt  shape.  For  further  details  of  the  theory  of  line  broadening,  the  reader  is 
referred  to  the  work  of  Griem  and  other  standard  literature. 


For  the  calculation  of  line  absorption  coefficients,  the  oscillator  strengths  for  all  transitions  have 
to  be  known.  The  oscillator  strength  for  a  transition  i  -» i'  is 


fi^  = 
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In  all  cases  considered  here,  a  single  outer  electron  participates  in  the  transition.  Assuming  an 
LS  coupling  scheme,  the  line  strength  maybe  written  as  [21] 
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Here  the  first  bracket  [.]  is  shorthand  for 

[J,  J',L,L']  =  (27  + 1)(27'  +  1)(2Z  + 1)(2 L'  + 1)  , 

and  the  {.}  brackets  refer  to  Wigner  6-j  symbols  [21],  whose  values  can  be  evaluated 
analytically.  The  designations  /  (  /' )  refer  to  the  orbital  angular  momentum  of  the  outer  electron, 
with  />  being  the  greater  of  these;  L  (L')  denotes  the  total  orbital  angular  momentum,  Lc  being 
that  of  the  core  electrons;  and  J  (J')  denotes  the  total  angular  momentum  of  the  state  i  ( i'  ). 
Finally,  the  square  of  the  radial  integral  (/?"/ )  is  calculated  in  the  Coulomb  approximation  of 
Bates  and  Damgaard  [22]. 

For  isolated  lines  in  non-hydrogenic  systems,  the  broadening  by  ions  is  of  secondary  importance. 
Therefore,  broadening  by  ions  is  not  taken  into  account  here,  with  the  following  exception:  In 
the  case  of  neutral  atoms  of  carbon,  nitrogen,  and  oxygen,  Seaton  [23]  presented  a  fit  formula 
that  yields  results  within  a  factor  of  2  of  Griem’s  [19]  published  results.  This  formula  is  based 
on  the  previously  described  line  strengths  and  is  used  for  Cl,  NI,  and  01.  Stehle  [24]  has 
published  tables  for  Stark  broadened  hydrogen  lines.  Her  calculation  is  based  on  the  Model 
Microfield  Method  (MMM)  for  both  the  electronic  and  ionic  broadenings  and  gives  an  accurate 
description  of  the  line  profile  from  the  center  to  the  wings.  These  data  for  the  hydrogen  Lyman, 
Balmer,  and  Paschen  series  are  used  up  to  electron  number  densities  of  1019  cm"3. 

The  photoionization  cross  sections  are  calculated  in  the  classical  theory  by  Kramers.  This  theory 
was  developed  for  hydrogen  and  hydrogen-like  species,  which  means  that  for  the  low-ionized 
plasmas  in  ETC  applications  the  corresponding  equations  are  only  correct  for  hydrogen  itself. 
Better  agreement  for  multi-electron  systems  is  obtained,  especially  for  high  exited  states  (which 
can  be  assumed  to  be  hydrogen-like),  by  use  of  effective  quantum  numbers  or  the  effective 


5 


nuclear  charges,  respectively.  For  the  photoionization  cross  sections  from  the  ground  states  of 
carbon,  nitrogen,  and  oxygen  published  data  of  more  detailed  calculations  have  been  adopted.  In 
addition,  so  called  Biberman  or  Bibennan-Norman  factors,  which  describe  the  difference  of 
cross  sections  for  multi-electron  species  to  the  hydrogen-like  one,  have  been  used  for  carbon, 
nitrogen,  oxygen,  and  their  first  ionization  stages.  As  will  be  shown  in  section  3,  calculated 
cross  sections  are  in  good  agreement  to  the  published  data,  e.g.,  to  NI  cross  sections  given  by 
Wilson  and  Nicolet  [25]. 

2.3  Method  of  Partial  Characteristics 

The  basis  for  quantitative  studies  of  radiative  energy  transport  is  the  radiation  transport  equation 
(RTE).  Exact  solution  of  this  equation  takes  many  hours  of  CPU  time  because  thousands  of 
absorption  coefficients  for  the  continuum  and  the  spectral  lines  have  to  be  considered  for  every 
iteration  step  when  solving  the  energy  balance  equation.  Simple  tabulation  of  the  absorption 
coefficients  in  the  presence  of  spectral  lines  is  also  virtually  impossible.  Therefore,  approximate 
but  effective  methods  of  integration  of  radiation  characteristics  are  necessary.  One  integral 
method  that  gives  values  very  close  to  the  exact  solution  is  the  Method  of  Partial  Characteristics 
(MPC)  presented  by  Sevastyanenko  [26,  27].  In  this  method,  the  complete  description  of  the 
radiation  field  is  possible  by  use  of  precalculated  partial  intensities 

exp 

o 

Here,  A /(£  X)  is  the  contribution  to  the  intensity  at  point  X  at  temperature  Tx  and  pressure  px, 
emitted  at  point  £  with  its  local  conditions  77  and  p.  and  taken  into  account  absorption  along  the 
path  OX . 

The  radiation  intensity  at  point  X  then  easily  can  be  estimated  by  a  simple  integration  along  a  ray 
OX: 

X 

t(x)~  Ja/(2',X>4U  (13) 

0 

and  the  radiation  flux  at  X  is  obtained  by  integration  over  all  directions 

s(x)=  Jj7(x)Q7Q.  (14) 

4  71 

Using  the  MPC,  it  is  possible  to  generate  tables  of  radiation  quantities  for  rapid  calculations  of 
radiation  transfer  in  gas-dynamic  simulations  [5,  28].  A  more  detailed  summary  of  the  MPC  can 
be  found  in  an  earlier  paper  [6]  or  in  the  original  work  of  Sevastyanenko  [26,  27], 
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2.4  Three-Dimensional  (3-D)  Radiation  Transfer 

The  ARL  NSRG2  CFD  code  [15]  solves  the  two-dimensional  (2-D)/axisymmetric  and  time- 
dependent  Navier-Stokes  equations  by  use  of  a  finite -volume  method.  That  means  for  a  given 
grid  cell,  it  balances  the  transfer  of  mass,  momentum,  and  energy  from  and  to  all  neighbor  grid 
cells.  For  the  radiation  field,  we  have  to  know  not  only  the  contribution  from  the  neighbor  grid 
cells,  but  also  the  contribution  from  all  cells  having  a  line-of-sight  to  the  cell  of  interest.  For  this 
reason,  a  quasi-3-D  integration  method  has  been  developed,  which  integrates  the  partial  radiation 
intensities  of  the  whole  3-D  volume  for  every  grid  point  of  the  2-D  CFD  grid,  i.e.,  all  grid  points 
4'(r',  <7>,  z')  of  the  3-D  volume  having  a  line-of-sight  to  a  grid  point  X  (r,  z)  of  the  2-D  CFD  grid 
are  taken  into  account.  The  necessary  conditions  at  any  point  of  the  3-D  volume,  e.g.,  the 
temperature,  can  easily  be  estimated  from  the  2-D  grid,  if  the  radial  and  longitudinal  coordinates 
r'  and  z'  of  the  point  C  are  known.  A  detailed  description  of  this  quasi-3-D  integration 
submodule  will  be  given  elsewhere  [7],  This  3-D  integration  submodule,  together  with  some 
routines  for  the  interpolation  of  radiation  characteristics  from  the  databases,  are  the  only  parts  of 
the  overall  radiation  code  which  will  be  integrated  in  the  CFD  code. 


3.  Validation 


The  final  ARL  plasma  flow  and  PPI  simulation  code  will  contain  highly  sophisticated 
submodules  like  the  radiation  transport  module,  a  module  for  solving  the  required  rate  equations, 
an  ablation,  and  ignition  module,  etc.  The  validation  of  the  coupled  modules  is  a  difficult  part, 
for  which  reliable  experimental  data  have  to  be  available,  too.  A  better  way  than  to  validate  the 
overall  program  in  a  single  step  is  to  validate  the  different  submodules  as  stand-alone  programs 
first.  In  the  case  of  the  radiation  part,  the  results  of  some  comparisons  to  published  theoretical 
and  experimental  data  will  be  shown. 

3.1  Atomic  Data 

About  150  calculated  oscillator  strengths  and  about  100  calculated  line  widths  have  been 
compared  to  the  extensive  compilation  of  the  U.S.  National  Institute  of  Standards  and 
Technology  (NIST)  [29]  and  the  standard  reference  tables  of  Griem  [17]  for  the  atoms  and  ions 
of  C,  N,  and  O  (for  hydrogen  the  theories  are  exact).  Almost  all  values  are  within  a  factor  of  2  to 
the  references,  whereas  most  of  them  deviate  only  by  a  maximum  of  20%.  For  example,  Tables  1 
and  2  show  some  results  for  the  oscillator  strengths  of  neutral  carbon  Cl  and  the  Stark 
broadening  of  single  ionized  carbon  CII. 


7 


Table  1.  Comparison  of  some  calculated  oscillator  strengths  of  neutral  carbon  Cl  to  data  of  Wiese  et  al.  [29]. 


Transition 
nl  ->  n'V 

Multiplat 

gi-gf 

E;  Ef 

(cm*1)  (cm’1) 

fi-r 

(this  work) 

fi-r 

(Wiese  et  al.  T291) 

2p-  3s 

3p  _  3po 

3-5 

16.40-60393.14 

6.615E-02 

5.88e-02 

1-3 

0.00-60352.63 

1.592E-01 

1.40e-01 

5-5 

43.40  -  60393.14 

1.191E-01 

1.04e-01 

3-3 

16.40-60352.63 

3.981E-02 

3.56e-02 

3-1 

16.40-60333.43 

5.316E-02 

4.71e-02 

5-3 

43.40-60352.63 

3.982E-02 

3.56e-02 

2p-3s 

lS-'P° 

1-3 

21648.01  -61981.82 

1.251E-01 

9.40e-02 

3s  -4p 

3p°  _  3p 

1-3 

60333.43-81325.76 

5.247E-03 

3.44e-03 

3-5 

60352.63-81343.99 

2.294E-03 

1.54e-03 

3-3 

60352.63-81325.76 

1.298E-03 

8.04e-04 

3-1 

60352.63  -  8131 1.01 

1.648E-03 

1.22e-03 

5-5 

60393.14-81343.99 

4.041E-03 

2.72e-03 

5-3 

60393.14-81325.76 

1.269E-03 

9.94e-04 

3s -3p 

3p°-3d 

3-5 

60352.63-69710.66 

3.757E-01 

3.99e-01 

1-3 

60333.43  -  69689.48 

5.004E-01 

5.34e-01 

5-7 

60393.14-69744.03 

4.215E-01 

4.42e-01 

3-3 

60352.63  -  69689.48 

1.250E-01 

1.29e-01 

5-5 

60393.14-69710.66 

7.502E-02 

7.60e-02 

5-3 

60393.14-69689.48 

4.990E-03 

4.95e-03 

3s  -4p 

¥pO  _  lp 

3-3 

61981.82-80562.85 

5.006E-03 

8.08e-03 

3d -6s 

3D  -  3P° 

3-5 

69689.48  -  86369.60 

1.886E-04 

1.93e-04 

5-5 

69710.66-86369.60 

1.700E-03 

1.74e-03 

3-3 

69689.48-86331.63 

2.829E-03 

2.89e-03 

3-1 

69689.48-86321.94 

3.768E-03 

3.85e-03 

7-5 

69744.03  -  86369.60 

6.820E-03 

6.93e-03 

5-3 

69710.66  -  86331.63 

5.100E-03 

5.21e-03 

Note:  Many  thousands  values  for  the  atoms  and  ions  of  hydrogen,  carbon,  nitrogen,  and  oxygen  are  taken 
into  account  in  the  radiation  transport  code. 


In  Figure  2,  a  comparison  of  a  calculated  continuum  spectrum  for  neutral  nitrogen  to  results  by 
Armstrong  et  al.  (published  by  Wilson  and  Nicolet  [25])  is  shown.  Good  agreement  is  obtained 
and,  whereas  the  spectrum  of  Armstrong  shows  the  sharp  edges  at  the  photoionization 
thresholds,  the  present  results  using  the  occupation  probability  formalism  are  more  physical. 

To  get  good  and  reliable  experimental  data  is  often  as  difficult  as  it  is  to  calculate  them.  Wiese 
et  al.  [10]  published  experimental  data  of  emissivities  of  the  hydrogen  Balmer  series,  which  have 
often  been  used  for  validation  of  theory  and  simulation.  In  Figure  3,  current  simulation  results 
are  compared  to  these  measured  data  at  two  different  temperatures  and  electron  number 
densities.  In  addition  to  free- free  Bremsstrahlung  and  bound-free  recombination  of  HI,  the  photo 
attachment  to  the  negative  hydrogen  is  included  in  the  simulation.  Both  theoretical  and 
experimental  emissivities  are  absolute,  i.e.,  no  scaling  has  been  applied. 
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Table  2.  Comparison  of  some  calculated  Stark  broadened  widths  (FWHM)  for  single  ionized  carbon  CII. 


Transition 
nl  — >  n  T 

Multiplat 

A, 

(A) 

T 

104 

(K) 

Ne 

1016 

(cm‘s) 

Ya 

(A) 

Yb 

(A) 

Yc 

(A) 

2p  -  3s 

2P°  -  2S 

858 

10 

10.0 

6.699e*  3 

6.889e*  3 

1.650e*  2 

20 

10.0 

4.737e*  3 

4.871e*  3 

1.300e*  2 

3s  -  3p 

2S  -  2P° 

6580 

10 

10.0 

8.324e*  i 

8.371e*  1 

1.192e+0 

20 

10.0 

5.886e*  i 

5.919e*  i 

1.038e+0 

3s  -  4p 

2S  -  2P° 

2174 

20 

10.0 

2.130e*  i 

2.267e*  l 

2.960e*  l 

3p  -  3d 

2P°  -  2d 

7234 

20 

10.0 

8.669e*  i 

1.068e+0 

1.198e+0 

3p  -  4s 

2P-2S 

3920 

20 

10.0 

8.604e*  i 

6.492e*  i 

9.600e*  i 

3p  -  4d 

2P°  -  2d 

2747 

20 

10.0 

1.153e+0 

1.392e+0 

9.360e*  i 

3d  -  4p 

2d  -  2P° 

5890 

20 

10.0 

1.667e+0 

1.896e+0 

2.180e+0 

3d  -  4f 

2d  -  2f° 

4267 

20 

10.0 

1.259e+0 

2.571e+0 

1.690e+0 

3d  -  5p 

2d  -  2P° 

3361 

20 

10.0 

1.999e+0 

2.754e+0 

1.836e+0 

3d  -  5f 

2d  -  2f° 

2993 

10 

10.0 

3.567e+0 

4.919e+0 

4.560e+0 

20 

10.0 

3.060e+0 

5.24  le+0 

3.840e+0 

Comparison  made  by  method  described  in  text. 

'’Comparison  made  by  a  simpler  approximation  for  the  matrix  elements. 

Comparison  made  with  values  given  by  Griem  [17]  (core  configuration  is  2s2(*S)). 

Note:  Many  thousands  widths  for  the  lines  of  atoms  and  ions  of  hydrogen,  carbon,  nitrogen,  and  oxygen 
are  taken  into  account  in  the  radiation  transport  code. 


Figure  2.  Comparison  of  calculated  total  effective  continuum  cross  section  for  N1  at  15,000  K.  The 
results  of  Armstrong  et  al.  [30]  have  been  published  by  Wilson  and  Nicolet  [25]. 
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Figure  3.  Emissivities  of  hydrogen  (a)  Ne=  2.9  x  1016  cm’*  and  T  =  10,850K  and  (b)  Ne=  9.3  x  1016  cm’5 
and  T  =  13,140  K,  compared  to  experimental  data  of  Wiese  et  al.  [10].  The  thin  lines  are  the 
contribution  due  to  lines  and  the  continuum.  The  smooth  dissolution  of  the  line  into  the 
continuum  agrees  very  well  with  the  measurements. 


3.2  Net  Emission  Coefficients 

To  validate  the  integrated  radiation  characteristics  stored  as  databases,  net  emission  coefficients 
have  been  calculated  and  compared  to  published  results.  Lowke  [31]  has  shown  that  the 
averaged  radiation  intensity  of  an  isothennal  cylinder  with  radius  R  is  approximately  the  same 
as  for  an  isothermal  sphere.  By  using  the  simpler  formula  for  the  sphere,  the  net  emission 
coefficient  is  given  as 


%  =  \lyk'vexp(-k’vR)dv . 
o 

From  equations  12  and  13  it  can  easily  be  seen  that  the  net  emission  coefficient  can  be  estimated 
from  the  radiation  characteristics  database  if  we  use  an  isothermal  temperature  profile  with  Tx=  7T 

and  radius  X  =  R.  In  Figure  4,  a  comparison  of  net  emission  coefficients  for  high  temperature  air 

* 

and  SF6  is  shown. 


The  calculations  of  SF6  have  only  been  performed  for  validation  purposes;  there  is  no  interest  in  this  gas  for  ETC  applications. 
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Figure  4.  Net  emission  coefficients  of  (a)  an  isothermal  air  plasma  at  different  radii  and  (b)  an  isothermal 
SF6  plasma  of  radius  0.1  cm,  compared  to  results  of  Sevastyanenko  and  Soloukhin  [32], 
Gleizes  et  al.  [33],  Liebermann  and  Lowke  [34],  and  Aubrecht  and  Gross  [35].  All  calculations 
have  been  performed  for  atmospheric  pressure. 

Current  efforts  are  the  validation  of  the  3-D  model.  First  results  indicate  that  it  is  working  well, 
so  in  the  next  step  of  the  STAS  program  an  interface  between  the  RT  code  and  the  CFD  code 
will  be  set  up  and  tested. 


4.  Radiation  Transport  in  CFD  Simulations 


There  are  two  aspects  of  radiation  in  ETC  CFD  simulation:  first,  the  energy  transfer  by  radiation 
in  the  hot  plasma  and  gas,  and  second,  the  interaction  of  the  radiation  energy  that  is  transmitted 
to  the  capillary  wall  or  propellant.  The  current  ARL  CFD  codes  do  not  include  the  ohmic 
heating  in  the  energy  balance,  i.e.,  the  equation  that  is  solved  may  be  written  as  [15] 

^^  +  V-(pve)+p-Vv  +Vq  +  VS  =  0,  (15) 

dt 

with  the  density  p,  the  velocity  v  ,  the  energy  e,  the  pressure  p,  and  the  thennal  heat  flux  q  .  The 
divergence  of  the  radiation  flux  S  will  be  calculated  directly  by  use  of  the  MPC  without 
numerical  differentiation.  From  the  inner  hot  core  of  the  plasma  to  the  outer  cold  gas  zone,  this 
VS  term  will  change  its  sign:  whereas  the  radiation  is  mainly  responsible  for  energy  losses  in 
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the  hot  core,  a  lot  of  radiation  energy  will  be  reabsorbed  at  the  cold  boundary  and  the 
surrounding  gas,  leading  to  a  temperature  increase  in  this  regions.  Regarding  the  interface,  only 
the  3-D  integration  submodule  together  with  some  interpolation  routines  and  the  MPC  databases 
will  be  coupled  to  the  CFD  code.  In  every  iteration  and  time  step  when  the  energy  balance  is 
solved,  the  CFD  code  will  call  the  radiation  submodule  with  the  actual  pressure  and  temperature 
profile  (for  all  grid  points),  and  the  divergence  of  the  radiation  flux  will  be  given  back. 

When  taking  into  account  material  ablation  in  the  capillary  and/or  ignition  of  a  propellant,  in 
addition  to  the  VS  values  the  radiation  energy  flux  at  the  wall  (or  grains)  is  required.  Calling  the 
RT  code  for  a  second  time,  which  in  that  case  will  use  a  different  integration  algorithm  [6],  this 
can  be  done.  The  resulting  radiation  energy  will  then  be  used  for  PPI  studies,  in  new  pyrolysis 
laws  and/or  ablation  modules  in  the  CFD  simulation.  Regarding  propellants,  which  may  be 
transparent  for  parts  of  the  electromagnetic  spectrum,  one  might  be  interested  in  the  energy 
transfer  within  different  frequency  bands.  Using  the  presented  radiation  transport  code,  either 
the  total  energy  transfer  from  the  vacuum  ultraviolet  (VUV)  to  the  infrared  (IR)  or  by  radiation 
of  some  frequency  bands  of  special  interest  can  be  calculated. 


5.  Summary 


A  new  radiation  transport  code  has  been  described,  which  will  be  used  in  the  ARL  NSRG2  code 
for  the  plasma-propellant  interaction  simulation.  With  the  implemented  method  of  partial 
characteristics  it  is  possible  to  precalculate  databases,  which  allow  one  to  calculate  the  energy 
transfer  by  radiation  iteratively  at  the  same  time  when  solving  the  balance  equations  of  mass, 
momentum,  and  energy  in  any  CFD  code  in  an  acceptable  CPU  time.  Some  efforts  for  the 
validation  of  the  program  have  been  presented,  and  it  has  been  shown  that  the  current  version  of 
the  code  yields  results,  which  are  in  very  good  agreement  to  published  theoretical  as  well 
experimental  data.  In  addition,  a  3-D  integration  method  has  been  implemented,  which  allows 
one  to  use  the  code  for  calculations  in  any  cylinder-symmetrical  volume  (the  precalculated 
partial  intensities  stored  in  the  databases  are  independent  of  the  geometry);  the  validation  of  this 
part  is  almost  finished.  In  a  next  step,  absorption  of  the  most  important  molecules  will  be  added 
to  the  code  and  current  efforts  are  the  coupling  of  the  RT  code  to  the  ARL  NSRG2  code.  The 
focus  of  the  current  activities  is  on  plasmas  containing  molecules,  atoms,  and  ions  of  hydrogen, 
carbon,  nitrogen,  and  oxygen,  but  the  program  is  in  no  way  limited  to  this.  Finally,  the 
performance  of  some  laboratory  experiments  is  underway  to  validate  the  CFD  code  with 
incorporated  radiation  simulation. 

Until  now,  the  radiation  transport  code  consists  of  more  than  40  subroutines  and  6  input 
databases/files  for  the  energy  level,  quantum  numbers,  etc.  Finally,  a  graphical  user  interface 
(GUI)  has  been  developed,  which  allows  one  to  use  the  code  without  extensive  knowledge  of 
quantum  mechanics  and  plasma  physics. 
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